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摘 要 : 空 时 信号 的 在 线 重 构 问 题 可 归结 为 对 差分 平滑 的 时 变 图 信号 的 恢复 问题 。 对 于 该 避 优 化 问题 , 现 有 的 基于 梯度 下 
降 法 的 分 布 式 重 构 算法 在 优化 问题 的 海 森 矩 阵 条 件数 较 大 时 收敛 速度 极 慢 , 在 单个 观测 区 间 内 算法 最 大 迭代 次 数 受 限时 


重 构 误 差 较 大 。 针 对 该 问题 ,提出 了 一 种 基于 近似 牛顿 法 的 分 布 式 在 线 重 构 算法 。 首 先 通 过 子 图 划分 将 原 优化 问题 分 解 
为 一 系列 子 图 上 的 局 部 优化 问题 ,并 求 出 该 局 部 问题 的 解 ,然后 对 子 图 间 的 局 部 解 作 融 合 平均 计算 ,得 到 近似 的 全 局 最 优 


解 ,再 依据 近似 解 与 实际 最 优 解 之 间 的 差距 ,证 明 以 此 方式 求 得 的 子 图 划分 与 融合 矩阵 具有 稀 玻 性 , 且 可 作为 原 优化 问题 


的 海 森 逆 近似 矩阵 ,最 后 将 该 近似 矩阵 蔡 换 至 经 典 的 牛顿 法 欠 代 公式 ,并 利用 该 近似 矩阵 的 结构 化 稀 玻 性 实现 分 布 式 运 
厂 。 念 真 结果 表明 ,与 现 有 算法 相 比 ,该 算法 收敛 速度 更 快 , 重 构 误差 更 小 ,所 需 通 信和 量 更 少 。 
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stract; The reconstruction problem of spatio-temporal signals can be cast as recovering differential smooth time-varying 
'graph signals. For the optimization problem, the existing distributed algorithm based on gradient descent method shows 


Slow convergence when the condition number of the Hessian matrix of the problem is large which leads to a large reconstruc- 


n error when the maximum iteration number is limited in an observation interval. Therefore. an online distributed recon- 
uction algorithm based on approximate Newton's method is proposed in the paper. whose principle is to decompose the 
ginal optimization problem into a series of local problems on subgraphs through subgraph decomposition and find these 


Olutions, and then obtain the approximate global optimal solution via fusion average of local solutions between each sub- 
'graph. According to 


he gap between the approximate solution and the actual one. it can be proved that the decompostion 


and fusion matrix obtained in this way is sparse and can be regarded as the approximate Hessian inverse. Hence, the algo- 


rithm replaces the approximate matrix into the classical Newton iterative formula which can be implemented in a distributed 


manner due to the structural sparsity of the approximate matrix. Simulation results show that the proposed algorithm has 


faster convergence rate and smaller reconstruction error. and requires less communication cost compared with the existing 


algorithm. 
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decomposition 


当今 信息 时 代 , 人 类 社会 无 处 不 充斥 着 海量 的 空 。 信息、 交通 网 络 中 的 车 流量 数据 和 生物 神经 元 网 络 中 
时 信号 。 随 着 科学 技术 的 进步 ,这 些 信号 逐渐 向 规模 ”的 生物 电信 号 等 。 现 实 世界 的 空 时 信号 因 能 量 受 


更 大 、 纬 度 更 高 、 结 构 更 复杂 的 趋势 发 展 , 如 无 线 传 感 限 、 环 境 污染 、 机 需 故 障 、 噪 声 干 扰 等 因素 影响 ,可 能 
器 网 络 中 的 温度 或 湿度 数据 、 社 交 媒 体 中 的 个 人 爱好 ” 是 缺失 或 不 可 靠 的。 以 无 线 传感器 网 络 中 的 监测 工 
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作为 例 , 由 于 电量 ,内存 和 处 理 能 力 有 限 , 传 感 器 节点 
无 法 做 到 对 数据 时 刻 观 测 , 因 此 整个 网 络 的 观测 数据 
可 能 是 不 完整 的 = 。 此 外 ,传感器 节点 还 可 能 受 硬 
件 老化 .噪声 干扰 或 剧烈 环境 变化 等 因素 影响 而 导致 
其 监测 数据 异常 呈 ,此 时 该 节点 的 可 靠 信 号 值 是 缺失 
的 。 在 空 时 信号 存在 数据 缺失 时 ,为 了 保证 后 续 信号 
处 理工 作 的 正常 进行 , 需 利 用 空 时 信号 的 关联 特性 对 
fri ob RU, 

图 信和 号 处 理 理 论 将 传统 信和 号 处 理 邻 域 中 的 传 里 
叶 变 换 、 滤 波 、 调 制 、 卷 积 等 概念 推广 到 非 规则 的 图 
上 ,是 研究 非 规则 域 信 号 的 有 力 工具 。 空 时 信号 
在 空间 域 与 时 间 域 上 的 关联 性 一 般 可 由 时 变 图 信号 
的 差分 平滑 性 来 描述 ""… 。 图 信号 的 差分 平滑 性 是 
指 该 信号 及 其 变化 量 在 图 上 相 邻 顶点 间 相 近 的 特性 。 
例如 ,在 温度 传 感 融 网 络 中 ,地 理 距 离 越 菲 近 的 节点 ， 


空 困 信号 进行 重 构 , 且 该 模型 具有 计算 复杂 度 低 、 重 
梅 醒 延 小 .可 分 布 式 求解 的 优点 5 。 然 而 , 现 有 的 求 
解 诬 模型 的 分 布 式 算法 均 基 于 梯度 下 降 法 po ,其 收 
伊 速 度 相对 缓慢 , 且 易 受 优化 问题 的 条 件数 影 
响 s 5 。 对 于 分 布 式 算法 ,过 慢 的 收敛 速度 会 导致 分 
布 或 系统 内 节点 通信 量 大 增 哺 ,从 而 缩短 系统 使 用 
寿 丛 ,因此 需要 设计 一 种 新 的 收敛 速度 更 快 且 更 稳定 
的 芬 布 式 算法 。 

泛 针 对 现 有 算法 "收敛 速度 相对 慢 的 问题 ,提出 
了 蝶 种 基于 近似 牛顿 法 的 空 时 信号 分 布 式 在 线 重 构 
算 活 。 通 过 子 图 划分 求 出 图 信号 重 构 优 化 问题 的 局 
部 解 ,再 进行 子 图 融合 ,得 到 近似 的 全 局 最 优 解 ,并 证 
明了 以 此 方式 求 得 的 子 图 划分 与 融合 矩阵 具有 稀 政 
性 ,可 作为 原 优化 问题 的 海 森 着 近似 ,最 后 将 近似 矩 
阵 代替 至 牛顿 法 迭代 公式 ,从 而 得 到 分 布 式 近似 牛 
顿 法 。 


1 空 时 信号 在 线 重 构 模 型 


1.1 图 信号 处 理 基 础 知识 


任意 空间 分 布 的 N 节点 网 络 可 建 模 为 一 个 图 。 
图 G 是 由 一 组 顶点 与 连接 顶点 的 边 构成 的 非 线 性 数 
Je 8H uy GS=V, E,W), 其 中 : V 为 顶点 集合 ; E 
为 边 集合 ; W c RVU 为 图 的 加 权 令 接 和 矩阵 。 一 般 
地 ,图 加 权 邻 接 矩 阵 可 完全 描述 一 个 图 的 拓扑 结构 ， 
其 权 值 满足 条 件 : 若 图 上 顶点 ; 与 顶点 7 相连 , 则 Wi > 
0; 否则 , W, —0, 对 于 相 邻 的 两 顶点 ! 和 7 , 顶点 间 
关联 性 或 相似 性 越 强 , 则 边 的 权 值 W; 越 大 。 除 了 加 


源 等 一 种 空 时 信号 的 分 布 式 在 线 重 构 算法 000000, 35 


CNMIMAaATiv 合 TFARHT! 
BUBBEARBVEP, xen] HR EBGEHWHTDBIEL-—D—W 
描述 图 的 拓扑 结构 ,其 中 D 为 图 的 度 和 矩阵 , 其 元 素 
满足 


N 
di = DW,. (D 


对 于 NN 顶点 的 图 G, 图 信号 可 定义 为 由 其 所 在 
图 G 的 顶点 集合 V 到 实数 集 R 的 映射 , 即 x :V 一 
R。 若 图 的 顶点 顺序 按 某 种 方式 确定 , 则 图 信号 可 简 
记 为 一 个 N 维 向 量 x — [x12 rs] ERY, H 
P r; 为 图 信号 在 顶点 i 上 的 信号 值 。 通 常 将 一 段 时 
间 工 内 随时 间 发 生变 化 的 图 信和 号称 为 时 变 图 信号 ， 
并 将 其 表示 为 由 一 组 各 时 刻 图 信号 向 量 组 成 的 矩阵 
和 一 [xiyx s xr] E RT, 
图 信号 的 平滑 性 是 指 图 信号 在 图 上 相 邻 两 顶点 
间 的 信号 值 差 异 较 小 的 特性 。 一 般 地 ,可 由 图 拉 普 拉 
斯 二 次 型 x'Lx 来 度量 图 信号 的 平滑 度 。 根 据 图 拉 
普 拉 斯 矩阵 的 定义 , 拉 普 拉 斯 二 次 型 可 改写 为 
rr 一 >) MW; —2;». (2) 


由 式 (2) 可 知 , 拉 普 拉 斯 二 次 型 的 值 越 小 , 则 权 值 越 大 
的 边 相 连 的 两 顶点 间 信 号 值 差异 越 小 ,因此 图 信号 越 
平滑 ;反之 , 则 图 信号 越 不 平滑 。 


1.2 差分 平滑 图 信号 在 线 重 构 模 型 


将 不 完整 的 空 时 信号 抽象 为 一 段 时 间 T. 内 经 过 
采样 的 时 变 图 信号 。 时 变 图 信号 的 采样 观测 模型 为 
y,—S,x,,d- n, r-—1.2.; T, (3) 
其 中 : y, 为 t+ 时 刻 图 信号 的 实际 观测 值 ; x*,。 为 图 信 
号 的 真实 值 ; n, 为 独立 同 分 布 的 观测 噪声 ; S, 为 采 
样 操作 。 采 样 操作 S, 是 一 个 对 角 和 矩阵 ,其 对 角 元 素 
T Ie EDU i 被 采样 , 则 其 第 i 个 对 角 元 素 为 5,,; = 
1; 否则 , s,,; 二 0。 经 过 采样 操作 ,图 信号 的 实际 观测 
值 存在 零 元 素 , 即 观测 数据 存在 缺失 的 情况 。 
对 于 数据 缺失 的 平滑 图 信号 ,可 通过 和 常用 的 平滑 
图 信号 在 线 重 构 模型 对 不 完整 信号 进行 恢复 : 


加 | 2 
min 5 | S,x, — y, ES (4) 


其 中 ,目标 函数 第 一 项 为 数据 匹配 项 ,人 迫使 重 构 信 和 号 
在 采样 点 靠近 真实 值 ,而 第 二 项 为 平滑 正则 项 ,用 来 
提升 重 构 信 号 的 平滑 度 。 虽 然 式 (4) 可 用 于 空 时 信号 
的 重 构 任 务 ,但 是 该 模型 只 考虑 了 图 信号 在 图 顶点 域 
上 的 平滑 性 ,而 未 充分 利用 信号 在 时 间 域 的 关联 性 ， 
因此 重 构 误差 较 大 "中 。 为 此 ,文献 [10] 指 出 , 空 时 信 
号 不 仅 本 和 映 在 空间 域 呈 平 滑 性 ,在 时 间 域 也 呈 平 滑 变 
化 ,并 提出 了 对 空 时 信号 重 构 效果 更 好 的 差分 平滑 图 
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言 号 在 线 重 构 横 型 : 


1 5.170 —- " 
in 2 | S,x, yı | 2 | 2 (x, X VULG =x "»P 


(5) 


Hep xa 为 上 一 1 时 刻 的 重 构图 信号 。 根 据 该 模型 ， 
当 已 知 当前 时 刻 观测 信号 和 上 一 时 刻 的 重 构图 信号 
时 ,可 利用 时 变 图 信号 的 差分 平滑 性 重 构 当前 时 刻 信 
号 值 。 当 t= 二 1 时, 令 ,1 二 0, 则 此 时 重 构 模型 (5) 
退化 为 模型 (4) ,可 见 模型 (5) 得 到 的 第 1 个 时 刻 重 构 
图 信号 是 平滑 的 , 且 所 有 时 刻 信 号 总 体 呈 差分 平滑 的 
特点 ,因此 该 在 线 重 构 模 型 充分 考虑 了 图 信号 时 间 顶 
点 联合 域 上 的 平滑 性 


2 LUNDI ENS 


一 虽然 文献 [10] 提 出 的 在 线 重 构 模型 (5) 可 对 空 时 
f 加 进行 有 效 重 构 ， 但 其 求解 算法 属于 一 阶 算法 , 收 
做 逮 度 相对 慢 ,上 且 易 受 优化 问题 海 森 矩 阵 的 条 件数 影 
咯 了 DD 对 此 ,设计 一 种 基于 近似 牛顿 法 的 求解 优化 模型 
(了 7 的 分 布 式 算法 ， 

入 模型 (5) 属 于 最 小 二 乘 优化 问题 ,可 分 别 求 出 其 
s g, GO TIR EABIE H, : 

k) —S,x, (k) — y, HOL GG) —3X 1); (6) 
=S, + 0L. (7) 
job k pu. 因为 图 拉 普 拉 斯 矩阵 工 是 一 个 
与 图 具有 相同 稀 玻 模式 的 局 部 矩阵 ， 所 以 梯度 g, Go) 
的 师 算 可 分 布 式 实现 ， 
三 mMk] n] (8) 
rr GO jJ; suubx, k]; — Lydi; + 
0 S, WjQv,GOd;—LDv, GO])0. (9) 


j€V,, 

其 中 , [Lg,(R)]; Juss g, CO 的 第 i 个 元 素 。V;, 为 

顶点 i 的 > 阶 邻 域 ,其 定义 为 :从 顶点 i 出 发 ， Er BE 

内 可 到 达 的 所 有 顶点 构成 的 集合 (包含 顶点 i 本 身 )。 

可 见 ,各 顶点 计算 梯度 信息 Lg, GO ]; 需要 与 其 一 阶 

邻 域 硕 点 交换 一 次 信息 Lv, C82 ]; . 
根据 凸 优化 理论 ,牛顿 法 属于 二 阶 算法 ,其 收 


1 


at AH 
iTGNfvV 呈 FS 月 T 


Re 


的 近似 矩阵 已 来 代替 海 森 

法 的 近似 迭代 过 程 : 

x, G +D =x, (k) t d, GO) —x, G) — Pg OD, 
(10) 


my 完成 牛顿 


其 中 d, Ck) 为 近似 牛顿 步 长 。 
2.1 EHE 


为 了 求 得 优化 问题 的 海 森 道 近 似 和 矩阵 P,, 首先 
将 图 G 划分 为 一 系列 重合 的 以 顶点 i 为 中 心 的 子 图 
AME 

G =U Gi; ， d1) 

Gin, = (Van oT a2) 

其 中 , G, 为 由 顶点 i 的 xi 阶 邻 域 顶 点 集合 及 

连接 这 些 顶 点 的 边 集合 EF;,, 组 成 的 子 图 。 根 据 此 子 

图 划分 方式 ,可 将 图 G 上 的 全 局 优化 问题 (5) 分 解 为 
一 系列 子 图 C,,， 上 的 子 优化 问题 : 


A 1 ri 
MIN g | SR, x, m" | 2 "p 


Re do: DUNS X Eu (13) 
优化 模型 (13) 中 , Rin 为 一 个 表示 局 部 操作 的 对 角 
矩阵 , 且 满 足 ; 若 图 G 上 的 顶点 7 ETK G: 内 , 则 
R:n 的 第 j 个 对 角 元 素 等 于 1; 否则 ,等 于 0。 依靠 
局 部 操作 R; > Pid 目标 函数 值 仅 与 子 图 
G.， 内 的 优化 变 , X, 相关 ,而 无 需 该 子 图 外 的 
其 他 变 tt, Dt. presen 一 个 局 部 优化 问 
题 ,可 根据 子 图 内 数据 求 出 局 部 最 优 解 : 

£,; = Rin HR D Ri, Qi 0L 1), (14) 
Hp CO! 表示 和 矩阵 求 伪 道 。 由 于 和 矩阵 R, HR. 
有 具 有 与 子 图 G;,, HRH. H R; HR 只 
有 对 应 于 子 图 Gn 项 点 的 行列 非 零 , 而 其 他 行列 均 
为 全 零 行 或 全 零 列 。 因 此 ,该 N 阶 和 矩阵 求 伪 逆 的 计 
算 过 程 可 在 子 图 内 局 部 进行 , 且 计 算 复杂 度 等 价 于 一 
D ri SERRÉ. 
直接 通过 式 (14) 求 得 的 局 部 解 损失 了 子 图 Gi 


外 的 信息 。 为 了 互相 弥补 各 子 图 内 信息 的 缺失 , 需 对 
子 图 的 局 部 解 作 信息 


敛 速度 较 快 , 且 对 海 森 矩阵 条 件数 不 敏感 ,因而 
可 满足 在 线 重 构 算法 对 收敛 速度 的 要 求 ,但 需 对 海 
森 和 矩阵 求 逆 。 和 矩阵 求 道 是 一 个 集中 式 运 算 过 程 , 因 
此 对 于 整个 图 而 言 ,无 法 分 布 式 计算 图 规模 大 小 的 
和 矩阵 逆 。 不 仅 如 此 ,矩阵 求 逆 的 计算 复杂 度 较 高 ， 
KHON’), 不 适用 于 大 规模 的 空 时 信号 重 构 任务 。 
为 了 保持 二 阶 算法 收敛 速度 的 优势 ,同时 又 可 实现 
算法 的 分 布 式 计算 , 需 寻 求 一 个 具有 结构 化 稀 玻 性 


EA! R2 EZ (15) 


EIL r2 |;év; " 

其 中 : e, 为 经 子 图 融合 后 得 到 的 对 原 优化 问题 (5) 的 
近似 解 ; | 了, | 为 集合 V; 的 非 零 元 素 个 数 ; x, 为 
控制 子 图 融合 规模 的 参数 ,通常 为 一 个 小 于 r 的 正 
整数 。 选 取 合适 的 参数 +,, 可 有 效 降低 子 图 融合 的 
边界 效应 59 。 
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综合 子 图 划分 (14) eee aid 
上 的 每 个 顶点 i 仪 需 先 与 其 x 阶 邻 域 顶 点 交换 一 
TB ORAS xu. HEB r 阶 邻 域 顶点 交换 一 Ra 
息 , 便 可 得 到 近似 最 优 解 , ,因此 整个 过 程 可 分 布 式 
实现 。 将 以 上 过 程 合并 写 为 整个 图 上 紧凑 的 和 矩阵 向 


量 乘积 形式 : 
t, OUR, bs pulo (16) 
结合 式 (14)， sr UI P, T 
P, = (ÙR; n) * D Ran Rin HR; Rs 
jEV | iEV 


(17) 

则 近似 解 e, 可 简写 为 

&, =P, (y, OLX 1). (18) 
因为 子 图 划分 与 融合 过 程 需 每 个 顶点 总 共 交 换 x 十 
; 芥 邻 域 顶 点 的 信息 ,所 以 矩阵 P, 的 第 i 行 仅 在 对 
MRR Gon, 顶点 的 位 置 上 存在 非 零 元 素 , 即 该 逢 
MRHAR BH 
CO 令 原 优化 问题 的 梯度 向 量 (6 为 零 ,可 得 原 问题 


e 


TTC : 


£^. I 


* =H (y, +0LX a. (19) 
Bam static. 近似 解闷 与 最 优 解 
的 BB 的 差异 源 于 近似 矩阵 P， "P 的 
id 


> l, — x ll; I PH,—I|; lx? las (20) 
sb I 为 单位 矩阵 。 根据 文献 [19- 20] , 随 着 参数 一 、 
ra Gb. | P,H, — I ||; E 
Ln 6 找到 一 对 足够 大 的 参数 ,使 PH, —Ils 
Tm 即 近似 解闷 , 与 最 优 解 x” 间 的 相对 误差 小 于 1. 
此 时 可 认为 矩阵 P, 是 海 森 逆 矩阵 HH 的 近似 。 


2.2. FARER 


选取 合适 的 参数 对 rm 、r* 后 ,利用 式 (17) 可 得 海 
森 道 近似 矩阵 也 ,， KEPUBPE AGRO $ 导 到 近似 牛 
顿 法 迭代 式 : 

x, (k +D) =x, k) (GR, x 
jEV 
)'R,., g, G2. (21) 
由 于 和 矩阵 P, RARD 3E X (D n p x OH, 
算法 的 分 布 式 计算 流程 如 下 : 

1) 输 入 数据 LX za. , Lyd: W; ,0; 

2) 初 始 化 [à (0) |]， =[X |]， 

循环 执行 : 

3) 计 算 [ v, (CR) |] ex. (CR) J; zz [ ; 


SIR.. Rin HR,., 
i€V 


页 域 顶 ， 点 交换 信息 D, "nob ,并 计算 
[g, C&) ]; za UE CR) ]; E Ly, |; 十 


4) 与 1 阶 领 


j€Via 
DI r 阶 领 域 顶 点 交换 信息 Dg, 00 T; 并 计算 
d, (k) ——[R;,, GS, - LR, TR; g, k); 


Oj r: Bra EATA [d,; Gk) J; 3f 


[ v, (&) J]; ); 


算 


ELO o 5 dak] | Vas Ta 
i€Vi,, 


DW De G +D]: Dr 003; + [4,0155 
直到 满足 迭代 终止 条 件 。 


8) 输 出 数据 Lx, li. 

在 分 布 式 近 似 牛 顿 算法 的 每 次 循环 中 ,各 节点 需 
先 与 1 阶 领域 顶点 通信 一 次 ,计算 梯度 Lg, (0 ];. $8 
后 与 7 阶 领域 项 点 通信 一 次 ,计算 局 部 牛顿 步 长 
d, (E), 55 r 阶 领 域 顶 点 通信 一 次 ,计算 近似 牛顿 
步 长 Lqd,(&)j;, 最 后 完成 近 代 。 大 量 仿真 结果 表明 ， 
在 空 时 信号 在 线 重 构 问 题 中 ,参数 选取 为 7 =r, —1 
时 ,该 近似 牛顿 算法 已 经 拥有 足够 快 的 收敛 速度 。 


3 仿真 结果 与 分 析 


通过 实验 对 比 本 算法 (参数 选取 为 71 =r =1 的 
分 布 式 近似 牛顿 算法 ) 与 文献 [10 算法 在 空 时 信号 在 
线 重 构 任务 中 的 性 能 。 仿 真 数 据 集 选用 人 工 合成 数 
据 集 和 海平 面 温度 数据 集 。 在 对 比 最 大 迭 代 次 数 开 
充足 情况 下 的 算法 性 能 时 ,性 能 评价 指标 设置 为 相对 
误差 : 


lx, GO — x; |， 
E, = : 22 
UO a ds um 
由 于 实际 采用 的 在 线 算法 通常 在 每 个 较 短 的 信号 观 


测 区 间 内 最 大 迭代 次 数 受 限 ,还 需 对 比 最 大 迭代 次 数 


K 受 限 情况 下 的 算法 性 能 ,此 时 性 能 评价 指标 设置 
为 累积 重 构 误 差 ; 
Ecr = ——X lx,— X, | 20 (23) 
2 REEEUE IIR CAR LE AR FEET ERE 
lx, k) —x,G — D ll; E 
lxi —DI, ec qt. (24) 
k >K. (25) 
实验 1 使 用 MATLAB 图 信号 处 理工 具 箱 


GSPbox! 沁 构建 一 个 100 顶点 的 随机 传 感 图 ,以 模拟 
现实 中 的 无 线 传感器 网 络 , 其 拓扑 结构 如 图 1 所 示 。 
根据 构建 出 的 图 随机 生成 一 组 100 个 顶点 .100 个 时 
刻 的 差分 平滑 时 变 图 信号 。 图 信号 生成 方式 为 : 先 随 
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机 生成 一 个 能 量 为 | xio l} 10! 的 低频 图 信号 , 作 
为 时 变 信号 在 第 1 个 时 刻 的 初始 信号 值 ,再 不 断 通 过 
Xano = Xio FL t 
生成 下 一 时 刻 的 图 信号 。 这 里 扰动 6, 是 一 个 零 均 值 
的 高 斯 白 噪 声 向 量 , 用 来 控制 信号 的 变化 量 , 以 保证 
合成 的 时 变 图 信号 的 随机 性 。 扰 动 6, 的 能 量 为 
6, 12 —1.3. 从 而 保证 了 合成 信号 满足 差分 平滑 
性 。 对 合成 的 图 信号 添加 噪声 ,并 进行 采样 ,以 模拟 
空 时 信和 号 实际 观测 数据 部 分 缺失 的 场景 。 图 信和 号 在 
各 顶点 的 观测 噪声 n, 设置 为 独立 同 分 布 的 高 斯 白 噪 
声 , 其 方差 为 0.01, 均 值 为 0。 采 样 方式 设置 为 在 各 
顶点 、 各 时 刻 独立 同 分布 的 随机 采样 ,采样 概率 为 
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1 随机 传 感 图 


Se E 2 为 最 大 和 迭代 次 数 充足 条 件 下 ,单一 时 刻本 算 
HO, =r, ==1) 与 文献 [10] 算 法 在 重 构 任务 中 相对 误 
差 写 算法 迭代 次 数 .通信 代价 的 关系 曲线 。 对 于 文献 
[363535 ,由 于 实验 中 优化 问题 的 海 森 矩阵 的 条 件数 
为 了 2. 3, 是 一 个 相对 小 的 数 ,文献 [10] 算 法 可 保持 一 
个 较 快 的 收敛 速度 ,但 依然 慢 于 本 算法 。 因 为 收敛 束 
度 较 慢 ,文献 [10] 算 法 达到 相同 误差 精度 所 需 通信 量 
也 比 本 算法 高 。 

考虑 迭代 次 数 受 限 的 情形 ,2 种 算法 在 整个 时 间 
段 内 的 累积 重 构 误 差 与 在 每 个 观测 间隔 [ — 1.21 内 
最 大 迭代 次 数 、 最 大 通信 和 量 之 间 的 关系 如 图 3 所 示 。 
由 图 3 可 见 , 当 最 大 迭代 次 数 较 小 ,如 开 — 1 时 ,文献 
[10] 算 法 由 于 收敛 速度 慢 ,无 法 在 每 个 时 刻 1 得 到 准 
确 的 重 构 信号 ,累积 误差 较 大 。 作 为 对 比 ,分 布 式 近 
似 牛 顿 算法 对 观测 间隔 内 算法 最 大 迭代 次 数 的 要 求 
并 不 高 ,即使 在 K =1 时 ,该 算法 的 误差 也 很 小 ,同时 
最 大 通信 量 需求 也 较 小 。 

实验 2 ”为 了 进一步 证 明 分 布 式 近似 牛顿 算法 
的 有 效 性 ,对 比 此 算法 Ci =r — D 与 文献 [10] 算 法 
对 太平 洋 海平 面 温度 数据 的 重 构 表现 。 数 据 集 选 取 


RT - - 一 文献 [10] 算 法 


和 迭代 次 数 
(a) 不 同 迭 代 次 数 下 的 相对 误差 


- -一 文献 [10] 算 法 
本 算法 


通信 代价 (x109) 
(b) 不同 通信 代价 下 的 相对 误差 


P 


图 2 最 大 迭代 次 数 充足 时 算法 性 能 对 比 


I2. 
1.0 一 一 一 文献 [10] 算 法 
AS 一 一 本 算法 
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最 大 迭代 次 数 
(a) 不 同 最 大 办 代 次 数 下 的 累积 重 构 误差 
2p 


1.0 L 一 一 一 文献 [10] 算 法 
Au —ÓO— 本 算法 
Z og 
z ; 
mx 0.6 
I 
B 
04 t 
0.2 : : ， 
0 0.5 1.0 1.5 2 
最 大 通信 代价 (x109 


(b) 不 同 最 大 通信 代价 下 的 累积 重 构 误 差 


图 3 最 大 迭代 次 数 受 限时 算法 性 能 对 比 


与 文献 L10] 相 同 ,为 太平 详 西 经 170" 至 90^. 南 纬 60° 
至 北纬 10° 传感器 网 络 温度 数据 ,其 拓扑 结构 如 图 4 
所 示 。 
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ri ^ 
2 P 
一 一 一文 献 L[10] 算 法 
L5r " 


180 W 150 W 120 W 90 W 60 W 


图 4 海平 面 温度 网 络 


同样 地 ,在 最 大 和 迭代 次 数 充足 的 条 件 下 ,2 种 算 
泪 杏 重 构 任务 中 相对 误差 与 迭代 次 数 ,通信 代价 的 关 
Ae npe 5 所 示 。 实 验 中 海 森 矩阵 的 条 件数 较 大 ， 
28-9 776 ,此 时 文献 [10] 算 法 的 收敛 速度 变 得 极为 组 
慢 - 从 而 导致 达到 目标 重 构 精 度 的 通信 需求 大 增 ,而 
代 调 式 近似 牛顿 算法 依然 保持 较 快 的 收敛 速度 。 先 
伐 次 数 受 限 时 ,2 种 算法 对 比如 图 6 所 示 。 由 于 收敛 
Bx 
次 精 度 的 重 构 信 号 ,而 分 布 式 近似 牛顿 算法 依然 有 较 
M o MuR, 


-一 文献 [10] 算 法 | 
一 一 木 算法 


= 100 
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XS DOR 
(a) 不 同和 迭代 次 数 下 的 相对 误差 
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信 代 价 (x1094) 
(b) 不 同 通信 代价 下 的 相对 误差 
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图 5 最 大 迭代 次 数 充 足 时 算法 性 能 对 比 


一 一 一 本 算法 


累积 重 构 误差 
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(a) 不 同 最 大 迭代 次 数 下 的 累积 重 构 误差 


一 一 文献 [10] 算 法 
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(b) 不 同 最 大 通信 代价 下 的 累积 重 构 误 差 


图 6 最 大 和 迭代 次 数 受 限时 算法 性 能 对 比 
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针对 现 有 空 时 信号 分 布 式 在 线 重 构 算法 收敛 速 
度 不 稳定 , 且 相 对 缓慢 的 问题 ,提出 了 一 种 基于 近似 
牛顿 法 的 空 时 信号 分 布 式 在 线 重 构 算法 。 通 过 子 图 
划分 求 出 图 信号 重 构 优 化 问题 的 局 部 解 ,再 进行 子 图 
融合 计算 近似 最 优 解 ,从 而 得 到 一 个 具有 稀 玻 性 的 海 
和 森 闭 近似 矩阵 ,最 后 将 该 近似 矩阵 代 蔡 至 牛顿 法 迭代 
公式 ,推导 出 分 布 式 近 似 牛 顿 法 。 仿 真 结果 表明 ,分 
布 式 近似 牛顿 算法 收敛 速度 更 快 ,通信 和 需求 更 低 , 且 
不 易 受 优化 问题 条 件数 的 影响 ,因此 更 适用 于 空 时 信 
号 分 布 式 在 线 重 构 任 务 。 本 研究 主要 是 在 已 有 优化 
模型 的 基础 上 进行 算法 的 改进 ,后 续 将 考虑 设计 新 的 
空 时 信号 在 线 重 构 模 型 , 以 进一步 降低 信号 重 构 
误差 。 
参考 文献 : 

[1] ORTEGA A,FROSSARD P, KOVACEVIC J, et al. 
Graph signal processing: overview, challenges, and ap- 
plications[ J ]. Proceedings of the IEEE, 2018, 106(5): 
808-828. 

[2] SANDRYHAILA,ALIAKSEI, MOURA, et al. Big da- 


ta analysis with signal processing on graphs: represen- 


tation and processing of massive data sets with irregu- 


lar structure[ J]. IEEE Signal Processing Magazine, 


40 


[3] 


[4] 


202382.00106v1 = 


aA 


chiriaXiv 


2014,31(52:80-90. 
YICK J, MUKHERJEE, GHOSAL D, et al. Wireless 
sensor network survey[ J]. Computer Networks. 2008, 
52(12) : 2292-2330. 

LORENZO P D. BANELLI P. BARBAROSSA S. et al. 
Distributed adaptive learning of graph signals[ J]. IEEE 


Transactions on Signal Processing. 2017.65(16) : 4193- 
4208. 

杨 杰 . 基于 图 信号 处 理 的 WSN 数据 异常 检测 与 恢复 
[D]. 桂林 :桂林 电子 科技 大 学 ,2019:1-4. 

杨 立 山 . 图 信号 采样 与 重建 研究 [Dj]. 北京 :北京 邮 有 
学 ,2018:9-12. 

韩 墨 . 图 信号 的 采样 与 重 构 理论 研究 [D]. 哈尔滨 : 哈 尔 
滨 工 业 大 学 ,2017:1-4. 

SHUMAN D I. NARANG S K, FROSSARD P. et al. 


The emerging field of signal processing on graphs: ex- 


大 


n 


tending high-dimensional data analysis to networks and 
other irregular domains[] ]. IEEE Signal Processing 
Magazine. 2013,30 (3) 183-98. 

SANDRYHAILA A, MOURA J. Discrete signal pro- 
cessing on graphs [J]. IEEE Transactions on Signal 
Processing, 2013,61(7) 11644-1656. 

QIU K, MAO X H, SHEN X Y, et al. Time-varying 
graph signal reconstruction[ J ]. IEEE Journal of Select- 
ed Topics in Signal Processing, 2017.11(6) :870-883. 
MAO X H.QIU K.LI T J.et al. Spatio-temporal signal 
recovery based on low rank and differential smoothness 
[J]. IEEE Transactions on Signal Processing. 2018. 66 
(23) :6281-6296. 

BOYD S, VANDENBERGHE L. Convex optimization 
[M]. Cambridge: Cambridge University Press, 2004; 
466-509. 


桂林 电子 科技 大 


[13] 


[14] 


[15] 


[16] 


[17 ] 


[18] 


[19] 


[20] 


[21] 


2 
TTF 


报 . 2023 年 2 月 
IH IV O ^ 
JAKOVETIC D. XAVIER J. MOURA J M F. Fast dis- 


tributed gradient methods[] ]. IEEE Transactions on 
Automatic Control.2014.59(5) : 1131-1146. 
MOKHTARI A.LING Q, RIBEIRO A. Network New- 
ton distributed optimization methods[J]. IEEE Trans- 
actions on Signal Processing, 2017 ,65(1) : 146-161. 
EISEN M. MOKHTARI A.RIBEIRO A. Decentralized 
Quasi-Newton methods[ J]. IEEE Transactions on Sig- 
nal Processing. 2017 ,65(10) :2613-2628. 

NARANG S K,GADDE A.SANOU E. et al. Localized 
iterative methods for interpolation in graph structured 
data[ C ]//IEEE Global Conference on Signal and Infor- 
mation Processing. Piscataway, NJ: IEEE Press. 2013: 
491-494. 

NOCEDAL J. WRIGHT S J. Numerical Optimization 
[M]. New York: Springer Science and Business Media. 
2006 :41-47. 

JIANG J Z, CHENG C, SUN Q Y. Nonsubsampled 
graph filter banks: theory and distributed algorithms 
[J]. IEEE Transactions on Signal Processing, 2019. 67 
(152:3938-3953. 

JIANG J Z. TAY D B. Decentralised signal processing 
on graphs via matrix inverse approximation[ J ). Signal 
Processing. 2019.,165(22) : 292-302. 

CHENG C.JIANG Y C.SUN Q Y. Spatially distribu- 
ted sampling and reconstruction[ J |. Applied and Com- 
putational Harmonic Analysis. 2019.47(1) 1109-148. 
PERRAUDIN N, PARATTE J, SHUMAN D. et al. 
GSPBOX: 
[EB/OL ]. (2016-03-15) [ 2021-03-20 ]. https: //arxiv. 
org/ pdf/1408. 5781. pdf. 


a toolbox for signal processing on graphs 


编辑 : 张 所 滨 


